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Abstract. A sequential dynamical system, or SDS, consists of an undirected graph Y, a vertex- 
indexed list of local functions g^y , and a word uj over the vertex set, containing each vertex at least 
once, that describes the order in which these local functions are to be applied. In this article 
we investigate the special case where F is a circular graph with n vertices and all of the local 
functions are identical. The 256 possible local functions are known as Wolfram rules and the re- 
sulting sequential dynamical systems are calledfinite asynchronous elementary cellular automata, 
or ACAs, since they resemble classical elementary cellular automata, but with the important dis- 
tinction that the vertex functions are applied sequentially rather than in parallel. An ACA is said 
to be ijj-independent if the set of periodic states does not depend on the choice of lo, and our main 
result is that for all n > 3 exactly 104 of the 256 Wolfram rules give rise to an w-independent 
ACA. In 2005 Hansson, Mortveit and Reidys classified the 11 symmetric Wolfram rules with this 
property. In addition to reproving and extending this earlier result, our proofs of w-independence 
also provide significant insight into the dynamics of these systems. 



Our main result, as recorded in Theorem 2.2, is a complete classification of the Wolfram rules 
that for all n > 3 lead to an a; -independent finite asynchronous elementary cellular automaton, 
or ACA. The structure of the article is relatively straightforward. The first two sections briefly 
describe how an ACA can be viewed as either a special type of sequential dynamical system 
or as a modified version of a classical elementary cellular automaton. These two sections also 
contain the background definitions and notations needed to carefully state our main result. Next, 
we introduce several new notations for Wolfram rules in order to make certain patterns easier to 
discern, and we significantly reduce the number of cases we need to consider by invoking the 
notion of dynamical equivalence. Sections 5 and 6 contain the heart of the proof. The former 
covers four large classes of rules whose members are cj-independent for similar reasons, and the 
latter finishes off three pairs of unusual cases that exhibit more intricate behavior requiring more 
delicate proofs. The final section contains remarks about directions for future research. 

1. Sequential Dynamical Systems 

Cellular automata, or CAs, are discrete dynamical systems that have been thoroughly studied 
by both professional and amateur mathematicians.' They are defined over regular grids of cells 
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'Stanislaw Ulam and John von Neumann were the first to study such systems, which they did while working at 
Los Alamos National Laboratory in the 1940s [9]. The German computer scientist Konrad Zuse proposed in 1969 
that the universe is essentially one big cellular automaton [13]. In the 1970s, John Conway invented the Game of 
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and each cell can take on one of a finite number of states. In addition, each cell has an update 
rule that takes its own state and the states of its neighbors as input, and at each discrete time step, 
the rules are applied and all vertex states are simultaneously updated. 

In the late 1990s, a group of scientists at Los Alamos invented a new type of discrete dy- 
namical system that they called sequential dynamical systems, or SDSs. In an SDS the regular 
grids used to define cellular automata are replaced by arbitrary undirected graphs, and the lo- 
cal functions are applied sequentially rather than in parallel. Their initial motivation was to 
develop a mathematical foundation for the analysis, simulation and implementation of various 
socio-technological systems [1, 2, 3, 4]. 

An SDS has three components: an undirected graph Y, a list of local functions ^Jy, and an 
update order uj. Start with a simple undirected graph Y with n vertices, label the vertices from 1 
to n, and recall that the neighbors of a vertex are those vertices connected to it by an edge. If F 
is a finite field and every vertex is assigned a value from F, then a global state of the system is 
described by an n-tuple y whose i"^ coordinate indicates the current state of the vertex i. The set 
of all possible states is the vector space F". 

Definition 1.1 (Local functions). A function F: F" ^ F" is called Y -local at i if (1) for each 
y G F", -F(y) only alters the z*^ coordinate of y and (2) the new value of the i^^ coordinate only 
depends on the coordinates of y corresponding to i and its neighbors in Y. Other names for such 
a function include local function and update rule. We use to denote a list of local functions 
that includes one for each vertex of Y. More precisely, = (Fi, F2, . . . , F„) where Fj is a 
function that is F-local at i. 

It is sometimes convenient to convert a local function F into another function with a severely 
restricted domain and range. 

Definition 1.2 (Restricted local functions). If z is a vertex with k neighbors in Y, then corre- 
sponding to each function F that is F-local at i, we define a function / : F^"*'^ F where the 
domain is restricted to the coordinates corresponding to i and its neighbors, and the output is the 
new value F would assign to the i*^^ coordinate under these conditions. It should be clear that F 
and / contain the same information but packaged in different ways. Each determines the other 
and both have their uses. Functions such as F can be readily composed, but functions such as / 
are easier to describe explicitly since irrelevant and redundant information has been eliminated. 

The local functions that are easiest to describe are those with extra symmetries. 

Definition 1.3 (Symmetric and quasi-symmetric rules). Let z be a vertex in Y with k neighbors, 
let F : F" F" be a F-local function at i and let / : F''^"^ ^ F be its restricted form. If the 
output of / only depends on the multiset of inputs and not their order, in other words, if the 
states of i and its neighbors can be arbitrarily permuted without changing the output of /, then 
/ (and F) are called symmetric local functions. If they satisfy the weaker condition that at least 
the states of the neighbors of i can be arbitrarily permuted without changing the output, then F 

Life, a two-dimensional CA, that was later popularized by Martin Gardner [5]. Beginning in 1983, Stephen Wolfram 
published a series of papers devoted to developing a theory of CAs and their role in science [7, 10, 11, 12]. This is 
also a central theme in Wolfram's 1280-page book A A^ew Kind of Science, published in 2002. 
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and / are quasi-symmetric. A list of local functions '^y is symmetric or quasi-symmetric when 
every function in the list has this property. 

The last component of an SDS is an update order. 

Definition 1.4 (Update orders). An update order a; is a finite sequence of numbers chosen from 
the set {1, . . . , n} such that every number 1 < i < n occurs at least once. If every number 
\ < i < n occurs exactly once, then the update order is simple. Let Wy denote the collection 
of all update orders and let Sy denote the subset of simple update orders. The subscript Y 
indicates that we are thinking of the numbers in these sequences as vertices in the graph Y . 
When considering an arbitrary update order, we tend to use the notation uj = {uji,uj2, . . . , ujm) 
with m = \uj\ (and m > n, of course), but when we restrict our attention to simple update orders, 
we switch to the notation tt = (tti, 71-2, • • • , 7r„). 

Definition 1.5 (Sequential dynamical systems). A sequential dynamical system, or SDS, is a 
triple {Y,^y,uj) consisting of an undirected graph Y, a list local functions ^y, and an up- 
date order u G Wy. If w is the sequence {lui,lu2, ■ ■ ■ , Um), then we construct the SDS map 
[dy,^] ■ IF" as the composition [g^y, u] := F^^ o ■ ■ ■ o F^^. 

With the usual abuse of notation, we sometimes let the SDS map [5^y , a;] stand in for the entire 
SDS. The goal is to study the behavior of the map [^y, u] under iteration. In this article we focus 
on the set of states in F" that are periodic and we use Per[5^y , u] C F*^ to denote this collection 
of periodic states. The set of periodic states is of interest both because it is the codomain of high 
iterates of the SDS map and the largest subset of states that are permuted by the map. 

Definition 1.6 (a; -independence). A list of F-local functions ^y is called uj -independent if 
Per[5'y,to'] = Per[^^y, cu'] for all update orders u.uo' G Wy and -k -independent if Per[^5^y,7r] = 
Per[5'y, tt'] for all simple update orders tt, tt' G Sy. 

Every cj-independent ^y is trivially vr-independent. More surprisingly, Reidys has shown that 
these two conditions are, in fact, equivalent. 

Theorem 1.7 ([8]). A list ^y of y-local functions is c<j-independent iff it is yr-independent. 

Even though -independence is too strong to expect generically, there are nonetheless many 
interesting classes of SDSs that have this property, including two classes where cj-independence 
is relatively easy to establish. 

Proposition 1.8. If for every simple update order vr G Sy, every state in Per[5'y, tt] is fixed by 
the SDS map \^y, tt], then ^y is yr-independent and thus a; -independent. 

Proof. If y is fixed by \^y, tt], then y is fixed by each Fi in ^y (the simplicity of tt means that 
were Fi to change the i^^' coordinate, there would not be an opportunity for it to change back). 
Being fixed by each Fi, y is also fixed by [5^y,a;] for all uj G Wy, which includes all of Sy. 
Since this argument is reversible, the SDS maps with simple update orders share a common set 
of fixed states. If, as hypothesized, these are the only periodic states for these maps, then ^y is 
TT-independent, and by Theorem 1 .7, a; -independent. □ 

In our second example, a; -independence is essentially immediate. 
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Proposition 1.9 (Bijective functions). If every local function Fj in is a bijection, then for 
every update order u E Wy, Per[^y, cj] = F". As a consequence ^y is a; -independent. 

Proof. Since every Fi is a bijection, so is the SDS map [^y, u] and a sufficiently high iterate is 
the identity permutation. □ 

This last result highlights the fact that cj-independence focuses on sets rather than cycles, since 
c<j-independent SDSs with different update orders quite often organize their common periodic 
states into different cycle configurations. In fact, the restrictions of cj-independent SDS maps 
with different update orders to their common periodic states can be used to construct a group 
encoding the possible dynamics over this set [6]. 

Collections of a; -independent SDSs also form a natural starting point for the study of stochas- 
tic sequential dynamical systems. Stochastic finite dynamical systems are often studied through 
Markov chains over their state space but in general this leads to Markov chains with exponen- 
tially many states as measured by the number of cells or vertices. For a; -independent SDSs one 
is typically able to reduce the number of states in such a Markov chain significantly, at least when 
focusing on their periodic behavior. 



2. Asynchronous Cellular Automata 

Some of the simplest (classical) cellular automata are the one-dimensional CAs known as 
elementary cellular automata. In an elementary CA, every vertex has precisely two neighbors, 
the only possible vertex states are or 1, and all local functions are identically defined. Since 
every vertex has two neighbors, the underlying graph is either a line or a circle and the restricted 
form of its common local function is a map / : F'^ F where F = F2 = {0, 1} is the field with 
two elements. There are 2*^ = 256 such functions, known as Wolfram rules, and thus 256 types 
of elementary cellular automata. Even in such a restrictive situation there are many interesting 
dynamical effects to be observed. The focus here is on the sequential dynamical systems that 
correspond to these classical elementary cellular automata. 

Let Y = CirCn denote a circular graph with n vertices labeled consecutively from 1 to n, 
and to avoid trivialities assume n > 3. (The sequential nature of the update rules in an SDS 
makes infinite graphs such as lines unsuitable in this context.) Since these are the only graphs 
considered in the remainder of the article, we replace notations such as Wy or Sy with Wn and 
Sn, etc. In Circ„ we view the vertex labels as residue classes mod n so that there is an edge 
connecting ito i + 1 for every i. 

Definition 2.1 (Wolfram rules). Let Fj : F" F" be a Circ„-local function at i and let /j : F^ ^ 
F be its restricted form. Since the neighbors of z are z — 1 and z + 1, it is standard to list these 
coordinates in ascending order in F^. Thus, a state y e F" corresponds to a triple j/j+i) 
in the domain of /j. Call this a local state configuration and keep in mind that all subscripts 
are viewed mod n. In order to completely specify the function Fi it is sufficient to list how the 
i^^ coordinate is updated for each of the 8 possible local state configurations. More specifically, 
let (ui-i, yi, Vi+i) denote a local state configuration and let Zi, yi+i) be the local state con- 
figuration after applying Fj. The local function Fj, henceforth referred to as a Wolfram rule, is 
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completely described by the following table. 



(2.1) 



Vi-iyiVi+i 


111 


110 


101 


100 


Oil 


010 


001 


000 


Zi 




as 


as 


a4 


as 


0.2 


ai 


ao 



More concisely, the 2^ = 256 possible Wolfram rules can be indexed by an 8-digit binary number 
a7a6asa4a3a2aiao, or by its decimal equivalent k = Yli=o ^-^^^ ■ There is thus one Wolfram rule 
k for each integer < A; < 255. For each such n, k and i let Wolf^^''^ denote the Circ„-local 
function : F" just defined, let wolf denote its restricted form /« : F^ F, and let 

2IJo[f^^') denote the list of local functions (Wolf f\ Wolff \ . . . , Wolf,^f ^). We say that Wolfram 
rule k is cu-independent whenever QUolf^'''' is a; -independent for all > 3. 

For each update order uj there is an SDS (Circ„, 2IJolf^^\ a;) that can be thought of as an 
elementary CA, but with the update functions applied asynchronously (and possibly more than 
once). For this reason, such systems are called asynchronous cellular automata or ACAs. We 
now state our main result. 

Theorem 2.2. There are exactly 104 Wolfram rules that are cu-independent. More precisely, 
2no[f(f') is a;-independent for all n > 3 iff A; G {0, 1, 4, 5, 8, 9, 12, 13, 28, 29, 32, 40, 51, 54, 57, 
60, 64, 65, 68, 69, 70, 71, 72, 73, 76, 77, 78, 79, 92, 93, 94, 95, 96, 99, 102, 105, 108, 109, 110, 
111, 124, 125, 126, 127, 128, 129, 132, 133, 136, 137, 140, 141, 147, 150, 152, 153, 156, 157, 
160, 164, 168, 172, 184, 188, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 204, 205, 
206, 207, 216, 218, 220, 221, 222, 223, 224, 226, 228, 230, 232, 234, 235, 236, 237, 238, 239, 
248, 249, 250, 251, 252, 253, 254, 255}. 

The main result of [6] states that precisely 11 of the 16 symmetric Wolfram rules are u- 
independent over Circ„ for all n > 3. Theorem 2.2 significantly extends this result, reproving 
it in the process. In addition to identifying a large class of a; -independent ACAs, the proof also 
provides further insight into the dynamics of these systems at both periodic and transient states 
and thus serves as a foundation for the future study of their stochastic properties. We conclude 
this section with two remarks about the role played by computer investigations of these systems. 

Remark 2.3 (Unlisted numbers). The "only if" portion of this theorem was established experi- 
mentally. For each 4 < < 9, for each < < 255, and for each simple update order tt G Sn, 
a computer program written by the first and third authors calculated the set Per[22Jo[f,f ^, tt]. For 
each of the 152 values of k not listed above, there were distinct simple update orders that led to 
distinct sets of periodic states, leaving the remaining 104 rules as the only ones with the potential 
to be a; -independent for all n > 3. Moreover, since a counterexample for one value of n leads 
to similar counterexamples for all multiples of n, these 104 rules are also the only ones that are 
eventually Lj-independent for all sufficiently large values of n. Because these brute-force calcu- 
lations are explicit yet tedious they have been omitted, but the interested reader should feel free 
to contact the third author for a copy of the software that performed the calculations. 

Remark 2.4 (Computational guidance). These early computer-aided investigations also had a 
major impact on the "if" portion of the proof. Once the computer results highlighted the 104 
rules that were cj-independent for small values of n, we identified patterns and clusters among 
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Figure 1. Grid notation for Wolfram rules 



the 104 rules, which led to conjectured lemmas, and eventually to proofs that our conjectures 
were correct. The computer calculations thus provided crucial data that both prompted ideas and 
tempered our search for intermediate results. 

3. Wolfram Rule Notations 

Patterns among the 104 numbers listed in Theorem 2.2 are difficult to discern because the 
conversion from binary to decimal obscures many structural details. In this section we intro- 
duce other ways to describe the Wolfram rules that makes their similarities and differences more 
immediately apparent. 

Definition 3.1 (Grid notation). For each binary number k = 0706050403020109 we arrange its 
digits in a grid. The 8 local state configurations can be viewed as the vertices of a 3-cube and 
we arrange them according to the conventional projection of a 3-cube into the plane. See the 
left-hand side of Figure 1 . Next, we can place the binary digits of k at these positions as shown 
in the center of Figure 1. The boxes have been added because the local state configurations come 
in pairs. When a local function is applied, the states of the neighbors of i are left unchanged, 
so that the resulting local state configuration is located in the same box. We call this the grid 
notation for k. The grid notation for Wolfram rule 29 = 0001 1101 is shown on the left-hand side 
of Figure 2. 

Because grid notation is sometimes cumbersome to work with we also define a very concise 4 
symbol tag for each Wolfram rule that respects the box structure of the grid. 

Definition 3.2 (Tags). When we look at the grid notation for a Wolfram rule, in each box we see 
a pair of numbers, 11, 00, 10, or 01, and we encode these configurations by the symbols 1,0,-, 



and X, respectively. In other words '1' = 



1 1 


, '0' = 





i ; 

•> 


1 


, and 'x' 


1 



The symbols are meant to indicate that when the states of the neighbors place us in this box, the 
local function updates the i^^ coordinate by converting it to a 1, converting it to a 0, leaving it 
unchanged, or always changing it. We label the symbols for the four boxes pi, p2, Ps and as 
shown on the right-hand side of Figure 1 and we define the tag of k to be the string PiPsP2Pi ■ 
The numbering and the order of the pj's has been chosen to match the binary representation as 
closely as possible, with the hope of easing conversions between binary and tag representations. 
The process of converting Wolfram rule 29 to its tag Ox-1 is illustrated in Figure 2. 



Definition 3.3 (Symmetric and asymmetric). The middle row of the grid contains the positions 
where the states of the neighbors are equal and the top and bottom rows contain the positions 
where the states of the neighbors are different. We call the middle row the symmetric portion of 
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Figure 2. Converting Wolfram rule 29 = 0001 1101 into its tag x - 1 . 



the grid and the top and bottom rows the asymmetric portion. In the tag representation, the begin- 
ning and end of a tag describes how the rule responds to a symmetric neighborhood configuration 
and the middle of a tag describes how it responds to an asymmetric neighborhood configuration. 
With this in mind we call p^pi the symmetric part of the tag k = P4P3P2P1 and we call psP2 its 
asymmetric part. 

Table 1 shows the 104 cj-independent Wolfram rules listed in Theorem 2.2 arranged according 
to the symmetric and asymmetric parts of their tags. The rows list all 16 possibilities for the 
symmetric part of the tag while the columns list only 10 of the 16 possibilities for the asymmetric 
part since only these 10 occur among the 104 rules. In addition, each row and column label has 
a decimal equivalent, listed next to the row and column headings, that add up to k. In this format 
the benefits of the tag representation should be clear. Far from being distributed haphazardly, the 
cj-independent rules appear clustered together in large blocks. Table 1 reveals a lot of structure, 
but some patterns remain slightly hidden due to the order in which the rows and columns are 
listed. For example, there is a 4-by-4 block of bijective rules obtained by restricting attention to 
the four rows that show up in the last column and the four columns that show up in the last row. 

Proposition 3.4 (Bijective rules). Wolfram rules 51, 54, 57, 60, 99, 102, 105, 108, 147, 150, 
153, 156, 195, 198, 201 and 204 are cu-independent. 

Proof. The 16 rules listed have tags where each pi is either - or x. These (and only these) 
Wolfram rules correspond to bijective local functions and by Proposition 1.9 the ACAs these 
rules define are cu-independent. □ 



4. Dynamical equivalence 

In this section we use the notion of dynamical equivalence to reduce the proof of Theorem 2.2 
to a more manageable size. Two sequential dynamical systems {Y, ^y,^) and {Y, , co') defined 
over the same graph Y are said to be dynamically equivalent if there is a bijection H : ^ F" 
between their states such that H o [5^y,c<j] = [^Jy,a;'] o H. The key fact about dynamically 
equivalent SDSs, which is also easy to show, is that H establishes a bijection between their 
periodic states. In particular, iJ(Per[5^y, u;]) = Per[^Jy, u;']. Thus, if 5y is an a;-independent 
SDS and for each u E Wy there exists an u' E Wy such that {Y,^y,uj) and {Y,^'y,uj') are 
dynamically equivalent using the same function H, then ^y is also u;-independent. 

Although there are 256 Wolfram rules, many give rise to dynamically equivalent ACAs. In 
particular, there are three relatively elementary ways to alter an ACA to produce another one 
that appears different on the surface, but which is easily seen to be dynamically equivalent to the 
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Table 1 . The 104 a; -independent Wolfram rules arranged by the symmetric and 
asymmetric parts of their tags. 



original. These are obtained by (1) renumbering the vertices in the opposite direction, (2) sys- 
tematically switching all Is to Os and Os to Is, or (3) doing both at once. We call these alterations 
reflection, inversion and reflection-inversion of the ACA, respectively. The term reflection high- 
lights the fact that this alteration makes it appear as though we picked up the circular graph and 
flipped it over. We begin by describing the effect renumbering has on individual local functions. 

Definition 4.1 (Renumbering). The renumbering of the vertices we have in mind is achieved 
by the map r: Circ„ Circ„ that sends vertex i to vertex n + 1 — i. For later use we extend 
this to a map r: Wn Wn on update orders by applying r to each entry in the sequence. 
More specifically, if u = (tui, u;2, • • • , ^m), then r{uj) = (r(a;i), r(u;2), . . . , r{ujm))- Finally, 
on the level of states we define a map i?: F" — * F" that sends sends y = (yi, ?/2, . . . , yn) to 
(?/ni • • • 5 ?/2i 2/1)5 and we note that R is an involution. 

Definition 4.2 (Reflected rules). If the vertices of Circ„ are renumbered, rule Wolf is applied, 
and then the renumbering is reversed, the net effect is the same as if a different Wolfram rule 
were applied to the vertex r{i). Let i be the number that represents this other Wolfram rule. The 
differences between k and i are best seen in grid notation. The renumbering not only changes 
the vertex at which the rule seems to be applied, but it also reverses the order in which the 3 
coordinates are listed in the restricted local form. Only the asymmetric local state configurations, 
i.e. the top and bottom rows of the grid, are altered by this change so that the grid for i looks like 
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a reflection of the grid for k across a horizontal line. We call i the reflection of k and we define 
a map refl : {0, . . . , 255} {0, . . . , 255} with refl(A;) = ^. On the level of tags, the only change 
is to switch order of p2 and p-^, so, for example £=0 1-x is the reflection of A;=0-lx. 

In short, when £ = refl(A;), R o Wolfj'^^ o R = Wolf|,^^^.-| and, since R is an involution, this can 
be rewritten as i? o Wolffs = Wolf o R. 

Proposition 4.3. If £ = refl (A;), then SUolf^'^-' is oj -independent iff SHolf^^^ is cj-independent. 

Proof. The value of ^ was defined so that i?o Wolf) = y^o\r^^-^oR. As a result, for any uj G Wn, 
the ACA (Circ„, Wol^^^\uj) is dynamically equivalent to the ACA (Circ„, Wol^^^\r{u)) since 

R o [2no [fi'^' , cu] = Ro Wolffs o • ■ ■ o Wolf (J) o Wolf^f) 
= Wolf^? ,o...oWolf^? ,oWolf7 

= [Woif^\r{oo)]oR. 

The argument at the beginning of the section now shows that the a; -independence of Wolfl^^ 
implies that of Woifl^\ but since i = ref\{k) implies k = refl(^), the converse also holds. □ 

Similar results hold for inversions as we now show. 

Definition 4.4 (Inverting). Let 1 and denote the special states (1, 1, . . . , 1) and (0, 0, . . . , 0) in 
F". Since the function i{a) = 1 — a changes 1 to and to 1, the map /: F" ^ F" sending y to 
1 — y, has this effect on each coordinate of y. The map / is an involution like R, and from their 
definitions it is easy to check that they commute with each other. 

Definition 4.5 (Inverted rules). If the states of Circ„ are inverted, rule Wolf,-^'' is applied, and then 
the inversion is reversed, the net effect is the same as if a different Wolfram rule were applied at 
vertex i. Let i be the number that represents this other Wolfram rule. The differences between 
k and £ are again best seen in grid notation. The pre-inversion of states effects the local state 
configurations as though the grid had been rotated 180°. The second inversion merely changes 
every entry so that Is becomes Os and Os become Is. Thus the grid for i can be obtained from 
the grid for k by rotating the grid and altering every entry. We call i the inversion of k and define 
a map inv: {0, . . . , 255} {0, . . . , 255} with \m(k) = L On the level of tags, there are two 
changes that take place. Boxes pi and p^ switch places as do boxes p2 and pa, but in process the 
boxes are turned over and the numbers changed. If we look at what this does to the entries in 
a box, 11 becomes 00, 00 becomes 11, while 10 and 01 are left unchanged. To formalize this, 
define a conjugation map c: {1, 0, — , x} {1, 0, — , x} with c(l) = 0, c(0) = 1, c(— ) = — , and 
c(x) = X. When k has tag P4P3P2P1, ^ has tag c{pi)c{p2)c{p'i)c{p4), so, for example, ^ = xO-1 is 
the inversion of A; = - 1 x. 

In short, when I = inv(A;), / o Wolf^-'^' o / = Wolffs and, since / is an involution, this can be 
rewritten as / o Wolff^^ = Wolffs o /. 

Proposition 4.6. If i = inv(A;), then Woifl^^ is cij -independent iff Wolf^^ is a; -independent. 



10 



M. MACAULEY 



J. MCCAMMOND 



H.S. MORTVEIT 





Fo 




















P2 
























8 










X 




1 o o 


132 


140 


TiA 71-1 




72 


24 


18 

J- u 


n 

u — 


A 


A 


1 




1 ^0 


ZU'f 


1 


1 


-0 


1 O O 

Izo 


1 O O 

12o 


136 


X- 


36 


1 AO 

108 


/^A 

60 


54 


n n 
u u 


n 
u 


n 
u 




o 


XX 


■^^ 

J J 


1 0^ 


J / 


CI 


1 
1 






141 


-0 


128 


200 


152 




01 


5 


5 


13 


10 


160 


232 


184 




-X 


129 


129 


137 


0- 


4 


76 


28 




Ox 


1 


1 


9 


01 


5 


77 


29 




1- 


164 


164 


172 


00 





72 






10 


160 


160 


168 


Ox 


1 


73 






xO 


32 


32 


40 





Figure 3. The 41 cj-independent Wolfram rules up to equivalence, separated 
into two tables by their behavior in asymmetric contexts. 



Proof. The value of (. was defined so that / o Wolf. = Wolf, o J. As in the proof of Proposi- 
tion 4.3 this implies that for any uj G Wn, the ACA (Circ„, WoV\\^\ u) is dynamically equivalent 
to the ACA (Circ„, 2nolf[f\ a;). The argument at the beginning of the section and the fact that 
£ = inv(A;) implies k = inv(^), complete the proof as before. □ 

As an immediate corollary of Propositions 4.3 and 4.6, when i = refl(inv(A;)) = \m{rei\(k)), 
22Jolf^'^^ is a; -independent iff 22Jo[f^^^ is a; -independent. If we partition the 256 Wolfram rules 
into equivalence classes of rules related by reflection, inversion or both, then there are 88 distinct 
equivalence classes and the 104 rules listed in Theorem 2.2 are the union of 41 of them. 

Figure 3 displays representatives of these 41 classes in pared down versions of Table 1. We 
used reflection and inversion to eliminate 5 of the 10 columns. Every rule with a 1 in the asym- 
metric portion of its tag is the inversion of a rule with a instead. In particular, the entries in 
the 3 columns headed -1,1- and 1 1 are inversions of the entries in the columns headed 0-, -0 
and , respectively. Next, since reflections switch p2 and p-^ we can also eliminate the columns 
headed -0, -x as redundant. This leaves the 5 columns headed 0, 0-, --, x- and xx. Since 
the last 3 do not contain Os or Is, further inversions, or inversion-reflections can be used to 
identify redundant rows in these columns. 

As mentioned above, the 41 rules listed in Figure 3 are representatives of the 41 distinct equiv- 
alence classes of rules whose u;-independence needs to be established in order to prove Theo- 
rem 2.2. The rows in each table have been arranged to correspond as closely as possible with the 
structure of the proof. For example, the first three rows of the table on the right-hand side are the 
9 equivalence classes shown to be a; -independent by Proposition 1.9. 
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Figure 4. Four major classes of a; -independent Wolfram rules. 
5. Major Classes 

In this section we prove that four large sets of Wolfram rules are oj-independent. All of the 
proofs are similar and, when combined with Proposition 1 .9, leave only 6 equivalence classes 
of Wolfram rules that need to be discussed separately. The main tool we use is the notion of a 
potential function. 

Definition 5.1 (Potential functions). Let F: X — > X be a map whose dynamics we wish to 
understand. A potential function for F is any map p: X — M such that p{F{x)) < p{x) for all 
X G X. A potential function narrows our search for periodic points since any element x with 
p{F{x)) < p{x) cannot be periodic: further applications of F can never return x to its original 
potential, hence the name. The only elements in X that are possibly periodic under F are those 
whose potential under p never drops at all. If we call the inverse image of a number in M a level 
set of p, then to find all periodic points of F, we only need to examine its behavior on each of 
these level sets. Finally, it should be clear that when non-decreasing functions are used in the 
definition instead of non-increasing ones, the effect is the same. 

Definition 5.2 (SDS potential functions). A potential function for an SDS such as (Y, 5^y, a;) 
is a map p: F" ^ M that is a potential function, in the sense defined above, for the SDS map 
[^Jy, u;]. The easiest way to create such a function is to find one that is a potential function for 
every local function Fj in Of course, p should be either a non-decreasing potential function 
for each Fi or a non-increasing potential function for each Fj, rather than a mixture of the two, 
for the inequalities to work out. When p has this stronger property we call it a potential function 
for since such a p is a potential function for (Y, ^y, for every choice of update order u. 

Proposition 5.3. Rules 0, 4, 8, 12, 72, 76, 128, 132, 136, 140 and 200 are a; -independent. 

Proof. If k is one of the numbers listed above, then its grid notation matches the leftmost form 
shown in Figure 4. (Each * is to be interpreted as either a or a 1 so that 16 rules share this 
form, the 1 1 listed in the statement and 5 that are equivalent to the listed rules or to previously 
known cases.) The 4 specified values mean that local functions never remove Os. Thus, the map p 
sending y G to the number of Os it contains is a non-decreasing potential function for 22Jo [fff ^ . 
Moreover, the local functions Wolf^'^^ cannot change y without raising p(y), so all periodic states 
are fixed states (for any update order), and by Proposition 1.8 Wolfl^^ is (jj-independent. □ 

For the next potential function, additional definitions are needed. 

Definition 5.4 (Blocks). A state y G F" is thought of as a cyclic binary n-bit string with the 
indices taken mod n, and a substring of y corresponds to a set of consecutive indices. We refer 
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to maximal substrings of all Os as 0-blocks and maximal substrings of all Is as 1-blocks. If 
a block contains only a single number it is isolated and if it contains more than one number 
it is non-isolated. The state y = 010110, for example contains one isolated 0-block and one 
non-isolated 0-block of length 2 that wraps across the end of the word. 

We study how these blocks evolve as the local functions are applied. The decomposition of a 
Wolfram rule into its symmetric and asymmetric parts is particularly well adapted to the study 
of these evolutions. The asymmetric rules either make no change or shrink a non-isolated 1- 
block or 0-block from the left or the right, depending on which of the 4 asymmetric rules we are 
considering. Similarly, the 4 symmetric rules either do nothing, they remove an isolated block or 
they create an isolated block in the interior of a long block. 

Proposition 5.5. Rules 160, 164, 168, 172 and 232 are cj-independent. 

Proof. If k is one of the numbers listed above, then its grid notation matches the second form 
shown in Figure 4. The specified values mean that (1) the only Os ever removed are the isolated 
Os and (2) isolated Os are never added. In particular, non-isolated blocks of Os persist indefinitely, 
they might grow but they never shrink or split, and the isolated Os, once removed, never return. 
Thus, the map p that sends y to the number of non-isolated Os in y minus the number of isolated 
Os in y is a non-decreasing potential function for Wolfl^\ As before, the local functions \No\ff^^ 
cannot change y without raising p(y), so all periodic states are fixed states (for any update order), 
and by Proposition 1.8 22Jo[f,^f ^ is cu-independent. □ 

Proposition 5.6. Rules 5, 13, 77, 133 and 141 are a; -independent. 

Proof. If k is one of the numbers listed above, then its grid notation matches the third form 
shown in Figure 4. This time the specified values mean that (1) the only Os that are removed 
create isolated Is, and (2) isolated Is are never removed and they never stop being isolated. Thus 
the map p that sends y to the number of Os in y plus twice the number of isolated Is in y is a 
non-decreasing potential function for Woifl^\ Once again, the local functions Wolf-'''' cannot 
change y without raising p(y), so all periodic states are fixed states (for any update order), and 
by Proposition 1.8 SUolfj-f^ is cj-independent. □ 

The argument for the fourth collection is slightly more complicated. 

Proposition 5.7. Rules 1, 9, 73, 129 and 137 are cj-independent. 

Proof. If k is one of the numbers listed above, then its grid notation matches the rightmost form 
shown in Figure 4. This time the specified values mean that (1) the only Os that are removed 
create isolated Is, but (2) isolated Is can also be removed. The map p that sends y to the number 
of Os in y plus the number of isolated Is in y is a non-decreasing potential function for 
but the difficulty is that there are local changes with p(Wolf •'^^(y) = p(y). This is true for the 
local change 000 010 and for the local change 010 000. All other local changes raise 
the potential, but the existence of these two equalities indicates that there might be (and there 
are) states that are periodic under the action of some SDS map [SHo tf^''-' , a;] without being fixed. 
Rather than appeal to a general theorem, we calculate its periodic states explicitly in this case. 
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Fix an update order oj G VF„ and, for convenience, let F: F" ^ F" denote the SDS map 
[2no[f,^f \ cj] : F" F". If ag = and y contains a substring of the form Oil, then p{F{y)) > 
p(y) and y is not periodic under F. This is because either (1) the substring remains unakered 
until its central coordinate is updated, at which point it changes to and p is raised, or (2) it is 
altered ahead of time by switching the 1 on the right to a (also raising p), or by switching the 
on the left to a 1 (impossible since ai = = 0). Analogous arguments show that if ag = and 
y contains the substring 110, or if ay = and y contains the substring 111, then y is not periodic 
under F. Let F be the subset of F" where these situations do not occur. More specifically, if 
03 = remove the states with Oil substrings, if ag = remove the states with 110 substrings, 
and if 07 = remove the states with 111 substrings. If all three are equal to 1, then P = F". 

We claim that P = Per[2iro[f[f\ u], independent of the choice of u. We have already shown 
P C Per[Wolfl^\u]. Note that P is invariant under F (in the sense that F{P) C P) since the 
allowed local changes are not able to create the forbidden substrings when they do not already 

exist. Moreover, F restricted to P agrees with rule 201 = x, the rule of this form with 

as = aQ = aj = 1, since whenever 03, ag or aj is 0, P has been suitably restricted to make this 
fact irrelevant. Finally, for every u rule 201 is bijective, thus F is injective on P, F permutes the 
states in P and a sufficiently high power of F is the identity, showing every state in P is periodic 
independent of our choice of u;. □ 

6. Exceptional Cases 

At this point there are only 6 remaining rules whose a; -independence needs to be established 
and they come in pairs: 28 and 29, 32 and 40, and 152 and 184. These final 6 rules exhibit 
more intricate dynamics and the proofs are, of necessity, more delicate. We treat them in order 
of difficulty. 

Proposition 6.1. Rules 32 and 40 are cj-independent. 

Proof. Let k be 32 or 40, let tt = (tti, 712, . . . , 7r„) G be a simple update order, and let 
p.^n denote the SDS map [SHolf^f^ vr] : F" F". The listed rules share the leftmost 

form shown in Figure 5 and it is easy to see that is the only fixed state (1 is not fixed and 
a2 = ttQ = means the rightmost 1 in any 1 -block converts to when updated). We also claim 
is the only periodic state of F. Once this is established, the cj -independence of SUolfj-f^ follows 
immediately from Proposition 1.8. 

The values ao = oi = 04 = mean non-isolated 0-blocks persist indefinitely, they do not 
shrink or split. Moreover, 02 = ag = means that each non-isolated 0-block adds at least one 
on its left-hand side with each application of F. In particular, any state y 7^ with a non-isolated 
0-block eventually becomes the fixed point 0. Thus no such y is periodic. 

The rest of the argument is by contradiction. Suppose that y is a periodic point of F other than 
and consider the z*'^ coordinates in y, F{y) and F(F(y)). We claim that at least one of these 
coordinates is and at least one of these is 1. This is because at least 4 out of the 5 local state 
configurations that do not involve non-isolated Os change the coordinate (and when k = 32 all 5 
of them make a change). The only way that yi does not change value in F(y) is if immediately 
prior to the application of Wolf,- , the local state configuration is Oil. Between this application 
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Figure 5 . Three final pairs of a; -independent Wolfram rules. 

of Wolfp^ and the next, the to the left is updated. It either is no longer isolated at this point 
(contradicting the periodicity of y) or it now becomes a 1. In the latter case, the application of 
Wolf •'^^ during the second iteration of F changes the z*^ coordinate from 1 to 0. Note that we 
used the simplicity of the update order to ensure that each coordinate is updated only once during 
each pass through F. Finally, suppose thati = tti and choosey, F(y) orF(F(y)) so the (i + l)^* 
coordinate is a 0. As soon as is applied, there is a non-isolated 0-block, contradicting the 

claim that y 7^ is a periodic point. □ 

Since it was easy to show that every state is periodic under the bijective Wolfram rule 156 (with 
tag -X--), we did not examine the evolution of its blocks. We do so now since its behavior is 
relevant to our study of the 4 remaining rules. 

Example 6.2 (Wolfram rule 156). Because the symmetric part of rule 156 is -- no isolated 
blocks are ever created or destroyed and thus the number of blocks is invariant under iteration. 
Moreover, the four values ai = 05 = and 02 = 03 = 1 mean that substrings of the form 01 
are fixed indefinitely, leaving the right end of every 0-block and the left end of every 1-block 
permanently unchanged. The other type of boundary can and does move since = x, and it is 
its behavior that we want to examine. Let tt G Sn be a simple update order and let F : F" ^ F" 
denote the SDS map [Wo{f^^^\Ti] : F". So long as y is not or 1, there is a 1-block 

followed by a 0-block and a corresponding substring of the form 01 ■ ■ ■ 10 ■ ■ ■ 01. (If y only 
contains one 0-block and one 1-block, then the first two digits are the same as the last two digits, 
but that is irrelevant here.) As remarked above, the beginning of the 1-block and the end of the 
0-block are fixed, but the boundary between them can vary. 

Suppose both blocks are non-isolated and consider the central substring 10 at positions i and 
i + 1. These are the only positions in the entire substring that can vary and the first one to 
be updated will change value. Assume the is updated first. The 1-block grows, the 0-block 
shrinks and the boundary shifts one step to the right. As we cycle through the local functions, the 
simplicity of tt guarantees that the [i + 2)'^^^ coordinate is updated before the (i + 1)*^* coordinate 
is updated a second time. Thus the boundary shifts one more step to the right. This argument 
continues to be applicable until the 0-block shrinks to an isolated 0. At this point, the is still 
updated before the 1 to its left is updated again, but this time the remains unchanged. When 
the 1 to its left is updated it changes back to a 0, the 1-block shrinks, the 0-block grows and the 
boundary shifts to the left. The same argument with left and right reversed shows that now the 
0-block continues to grow until the 1-block shrinks to an isolated 1, at which point the shifting 
stops and the boundary starts shifting back in the other direction. 

Proposition 6.3. Rules 152 and 184 are -independent. 
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Proof. Let A; be 152 or 184, let vr = (tti, 7r2, . . . , 7r„) G S'n be a simple update order, and let 
F : F" ^ denote the SDS map [Wo Ifi''^ , n] : F" F". The listed rules share the second form 
shown in Figure 5 and it is easy to see that and 1 are the only fixed states (since 02 = ^6 = 
means the rightmost 1 in any 1-block converts to when updated). We also claim and 1 are 
the only periodic states of F. Once this is established, the cj -independence of Wo{\^^^ follows 
immediately from Proposition 1.8. 

Since isolated blocks are never created, the map p that sends y to the number of blocks it 
contains is a non-increasing potential function for SHotf ^f-'. Moreover, since the only differences 
between rule 156 and rules 152 and 184 are that rule 152 removes isolated 1-blocks and rule 184 
removes both isolated 1-blocks and isolated 0-blocks, the map F agrees with [SUotfj^^^^'', tt] so 
long as it is not called upon to update an isolated 1-block (or an isolated 0-block when k = 184). 
The long-term behavior of rule 156, however, as described in Example 6.2, shows that under 
iteration every y not equal to or 1 eventually updates such an isolated block, removing it and 
decreasing p, thus showing that such a y is not periodic. □ 

Finally, the argument for Wolfram rules 28 and 29 is a combination of the difficulties found in 
the proofs of Propositions 5.7 and 6.3. 

Proposition 6.4. Rules 28 and 29 are -independent. 

Proof. Let k be 28 or 29, let vr = (tti, 7r2, . . . , vr^) G S'^ be a simple update order, and let 
F: F" ^ F" denote the SDS map [Wolfl^\ vr] : F" ^ F". The listed rules share the rightmost 
form shown in Figure 5 and the values 05 = and 02 = 1 mean that isolated blocks are never 
removed. Thus the map p that sends y to the number of blocks it contains is a non-decreasing 
potential function for Wolfl^\ The four values oi = 05 = and 02 = 03 = 1 mean that 
substrings of the form 01 persist indefinitely, as in Wolfram rule 156. In fact, so long as p is 
unchanged, the behavior of F under iteration is indistinguishable from iterations of the map 
[Wo[f'^^^\n]. Consider a substring of the form 01 ■ ■ ■ 10 ■ ■ ■ 01 and suppose that the length of 
the 1-block on the left plus the length of the 0-block on the right is at least 4. We claim that any 
y containing such a substring is not periodic under F. If it were, the evolution of this substring 
would oscillate as described in Example 6.2 and at the point where the 0-block shrinks to an 
isolated 0, the 1-block on the left contains the substring 111. Moreover, between the point when 
that penultimate becomes a 1 and the point when it is to switch back, the substring 111 is 
updated, increasing p. When k is 29, a similar increase in p can occur when the 1-block shrinks 
to an isolated 1 and the 0-block contains the substring 000. In neither case can a state containing 
a 1-block followed by a 0-block with combined length at least 4 be periodic under F. 

Next, note that when k = 29 both of the special states and 1 are not fixed, but that for k = 
28 1 is not fixed, while is fixed. Let P be the set of states containing both Os and Is that do 
not contain a 1-block followed by a 0-block with combined length at least 4, and, when k = 28, 
include the special state as well. Because we understand the way that such states y G P evolve 
under Wolfram rule 156 (Example 6.2), we know that at no point in the future does a descendent 
of y ever contain a substring of the form 111 or 000. Thus P has been restricted enough to make 
the values of aj and Oq irrelevant, and F sends P into itself. Moreover, since F agrees with 
[Snolf^^^^-* , tt] on P, and this map is injective, F is injective on P, F permutes the states in P and 
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Number of flips 


1 2 


3 4 5 


6 7 8 


Number of a; -independent rules 


1 8 26 


34 26 4 


4 1 


Number of rules 


1 8 28 


56 70 56 


28 8 1 


Percentage 


100% 100% 93% 


61% 37% 7% 


14% 0% 100% 



Table 2. The number of flips and the probability of cu-independence. 



a sufficiently high power of F is the identity, showing every state in P is periodic, independent 
of our choice of tt. Now that we know that Snotfj-f^ is 7r-independent, u;-independent follows 
from Theorem 1.7. □ 



7. Concluding Remarks 

Now that the proof of Theorem 2.2 is complete, we pause to make a few comments about it and 
the 104 cj-independent Wolfram rules it identifies. For each of the 8 local state configurations. 
Wolfram rule k either leaves the central coordinate unchanged or it "flips" its value. The number 
of local state configurations that are flipped in this way is strongly correlated with the probability 
that a given rule is (^^-independent. See Table 2. The numbers in the third row are the binomial 
coefficients (^), since they clearly count the number of Wolfram rules with exactly i flips. The 
key facts illustrated by Table 2 are that virtually all of the rules with at most 2 flips are uj- 
independent, the percentage drops off rapidly between 2 and 6 flips, and a;-independence is very 
rare among rules with 6 or more flips. In fact, all 5 such rules are a; -independent because they are 
bijective. It would interesting to know whether this observation can be quantitatively (or even 
qualitatively) extended to a rigorous assertion about more general SDSs. 

Next, there are two aspects of Theorem 2.2 that we found slightly surprising. First, we did not 
initially expect the set of rules that were cj-independent for small values of n to match exactly 
the set of rules that were a; -independent for all values of n > 3. The second surprise was that the 
during the course of the proof we found that the Wolfram rules truly are local rules, in the sense 
that their set of periodic points tended to have essentially local characterizations. 

Finally, although the focus of this article was solely the classification of the 104 cj-independent 
Wolfram rules, and not the dynamics of these rules per se, many interesting dynamical properties 
arose in the course of the proof. We are currently studying the dynamics and periodic sets for all 
256 Wolfram rules in greater detail, as well as examining how the sets of periodic states under 
an a; -independent Wolfram rule get permuted as the update order is altered. The latter situation 
involves an object called the dynamics group of an u; -independent SDSs. We plan on publishing 
these further results in a future article that builds on the results described here. 
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